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Abstract 

Proton  exchange  membrane  (PEM)  fuel  cells  have  been  under  development  for  many  years  and  appear  to  be  the  potential  solution  for  many 
electricity  supply  applications.  Modelling  and  computer  simulation  of  PEM  fuel  cells  have  been  equally-active  areas  of  work  as  a  means  of 
developing  better  understanding  of  cell  and  stack  operation,  facilitating  design  improvements  and  supporting  system  simulation  studies. 

In  general,  fuel  cell  models  must  be  able  to  predict  both  activation  and  concentration  polarizations  at  both  anode  and  cathode.  Normally  these 
predictions  require  values  of  the  concentration  of  the  reactant  gases  (i.e.  H2  and  02)  at  the  interface  between  the  catalyst  and  the  electrolyte. 
Electrolytes  of  interest  could  include  various  dilute  acids  or  polymeric  membranes  such  as  Nation™  so  that  gas  solubilities,  in  the  form  of  Henry’s 
Law  constants,  could  be  required  for  a  diverse  group  of  solvents. 

Published  solubility  data  have  been  evaluated  and  a  number  of  Henry’s  Law  correlations  are  proposed. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Over  the  years,  there  has  been  an  extensive  literature  on  the 
theory  and  practice  of  fuel  cells  in  general  and  proton  exchange 
membrane  (PEM)  fuel  cells  in  particular.  Several  good  broad 
references  exist  to  give  an  overview  of  this  literature  [1-5]. 
There  is  also  an  extensive  literature  on  modelling  of  PEM 
fuel  cells.  Much  of  this  has  been  cited  and  used  in  previous 
publications  from  our  group  [6-12]  and  was  recently  summa¬ 
rized  [13,14].  Our  modelling  work  has  primarily  dealt  with  the 
development  of  the  ‘Generalized  Steady  State  Electrochemical 
Model’  (GSSEM). 

In  general,  all  fuel  cell  models  must  be  able  to  predict  both 
activation  and  concentration  polarizations,  both  of  these  loss 
terms  generally  being  a  function  of  the  concentration  of  the  reac¬ 
tant  gas  in  the  electrolyte.  These  polarizations  would  generally 
require  prediction  at  both  the  anode  and  the  cathode.  For  a  PEM 
fuel  cell,  therefore,  these  predictions  would  normally  require 
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values  of  the  concentration  of  both  of  the  reactant  gases  (i.e.  H2 
and  O2)  at  the  interface  between  the  respective  catalyst  and  the 
electrolyte. 

For  the  initial  development  of  a  database  for  the  prediction  of 
cell  parameters,  research  results  are  typically  obtained  in  single 
electrode  studies  in  lab-type  electrochemical  cells.  Electrolytes 
of  interest  at  this  stage  could  include  various  dilute  acids  such  as 
H2SO4  or  HC1  and  possibly  concentrated  acids  such  as  H3PO4. 
The  subsequent  application  of  such  a  database  to  the  modelling 
of  a  PEM  fuel  cell  requires  knowledge  of  reactant  concentra¬ 
tions  in  polymeric  membrane-electrolytes  such  as  Nafion™. 
Gas  solubilities,  in  the  form  of  Henry’s  Law  constants,  could 
therefore  be  required  for  a  diverse  group  of  solvents. 

Published  solubility  data  have  been  evaluated  and  a  number 
of  Henry’s  Law  correlations  are  proposed.  Procedures  to  use 
such  correlations  are  summarized. 

2.  The  PEM  fuel  cell  (PEMFC) 

The  PEM  fuel  cell  has  been  extensively  described  in  the  lit¬ 
erature.  Briefly,  there  are  three  major  components  that  are  of 
interest  to  all  PEM  fuel  cell  models: 
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Nomenclature 

A 

Arrhenius  pre-exponential  in  Eq.  (6)  (A  cm-2) 

B 

Arrhenius  slope  in  Eq.  (6)  (K) 

c 

concentration  (mol  cm-3) 

cgas 

concentration  of  dissolved  reactant  gas  in  the  elec¬ 
trolyte  (mol  cm-3) 

Cm 

concentration  of  active  Pt  sites  in  sites  cm-2  or 
moles  of  sites  cm-2 

E 

thermodynamic  emf  (the  maximum,  equilibrium, 
voltage  theoretically  possible  for  the  particular 
temperature  and  reactant  partial  pressures) 

fe 

‘function  of  parameter  in  Eq.  (7)  (dimension¬ 

less) 

GSSEM 

Generalized  Steady  State  Electrochemical 
Model 

H 

Henry’s  Law  ‘constant’  (atm  cm3  mol-1) 

HOR 

hydrogen  oxidation  reaction  (Eq.  (1)) 

i 

current  density  (A  cm-2) 

io 

exchange  current  density  (A  cm-2) 

MEA 

membrane  electrode  assembly 

ORR 

oxygen  reduction  reaction  (Eq.  (2)) 

P 

gas  partial  pressure  (atm) 

PEM 

‘polymer  electrolyte  membrane’  or  ‘proton 
exchange  membrane’ 

PEMFC 

PEM  fuel  cell 

RH 

relative  humidity  (expressed  as  %) 

T 

temperature  (K) 

V 

voltage  that  appears  across  the  fuel  cell  terminals 
(V) 

Greek  letters 

0 

polarization  (i.e.  overvoltage  or  loss)  (V) 

0 

fractional  surface  coverage  (e.g.  of  chemisorbed 
hydrogen  atoms  for  the  HOR) 

Subscripts 

a 

anode 

av 

average 

act 

activation  or  actual 

c 

cathode 

ch 

channel 

cone 

concentration  (as  related  to  ‘concentration 
polarization’,  the  voltage  loss  due  to  slow  mass 
transfer) 

h2 

hydrogen 

H 

hydrogen  atom 

h3o+ 

hydrated  proton  (i.e.  H  +  H2O) 

i 

at  the  gas-aqueous  interface  or  species  (e.g.  H2, 
02,  H,  etc.) 

02 

oxygen 

ohmic 

related  to  ohmic  (iR)  losses 

o 

at  zero  polarization  and  zero  net  current  condition 
(i.e.  at  equilibrium) 

Superscripts 

*  at  the  interface  of  the  Pt  catalyst  and  the  surround¬ 

ing  aqueous  medium  (the  ‘reaction  interface’) 
sat  saturated  (e.g.  water  vapour  at  RH  of  100%) 


(i)  The  anode  is  the  electrode  where  the  fuel,  hydrogen,  is 
converted  to  protons  and  where  electrons  are  released  to 
travel  through  the  external  circuit  and  perform  electrical 
work.  The  overall  half-cell  reaction  is 

H2  =>■  2H+  +  2e~  (1) 

The  anode  reaction  occurs  in  the  presence  of  a  catalyst, 
typically  Pt  (although  mixtures  such  as  Pt/Ru  may  be  used 
if  there  are  carbon-containing  species  in  the  anode  feed). 
Dissolved  hydrogen  gas  must  be  present  adjacent  to  the 
catalyst  surface.  The  ‘hydrogen  oxidation  reaction’  at  the 
anode  is  commonly  referred  to  as  the  ‘HOR’. 

(ii)  The  membrane,  separating  the  anode  and  cathode,  is  a  thin 
layer  of  a  special  polymer,  commonly  Dupont  Nation™, 
a  perfluorinated  polymer  with  sidechains  terminating  in 
sulphonic  acid  groups.  The  H+  ions,  each  hydrated  with 
several  water  molecules,  pass  through  the  membrane  from 
the  anode  to  the  cathode,  with  the  membrane  serving  as  the 
“electrolyte”  of  the  electrochemical  cell. 

(iii)  The  cathode  is  the  electrode  where  the  oxidant,  oxygen 
from  a  pure-oxygen  or  from  an  oxygen-containing  stream, 
first  catalytically  recombines  with  the  electrons  returning 
from  the  external  circuit 

02  +4e“  202~  (2) 

and  then  permits  the  completion  of  the  overall  half-cell 
chemical  reaction 

2H+  +  O2"  H20  +  heat  (3) 

The  ‘oxidation-reduction  reaction’  at  the  cathode  is 
commonly  referred  to  as  the  ‘ORR’.  The  cathode,  typi¬ 
cally,  also  uses  Pt  as  the  catalyst.  Dissolved  oxygen  must 
be  adjacent  to  the  catalyst  surface. 

3.  The  essentials  of  a  fuel  cell  model 

The  shortest  form  of  a  general  fuel  cell  model  is  typically 
given  by: 

V  =  E  +  77 act, a  T  0 act, c  T  0 ohmic  (4) 

where  the  rj  terms  are  all  losses,  i.e.  they  are  negative  in  Eq. 
(4).  The  first  two  of  these,  77act,a  and  77act,c>  the  “activation  polar¬ 
izations”,  primarily  incorporate  the  chemical  kinetic  parameters 
that  determine  the  rates  of  reaction  at  the  anode  and  cathode, 
respectively.  The  prediction  of  these  activation  polarizations  typ¬ 
ically  involves  the  use  of  a  reaction  rate  parameter,  the  exchange 
current  density,  itself  a  function  of  the  concentration  of  reactant 
in  the  electrolyte.  The  third  loss  term,  T/ohmic,  deals  with  the  “iR” 
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electrical  loss  associated,  primarily,  with  the  flow  of  hydrated 
H+  ions  through  the  membrane  electrolyte. 

Depending  on  the  complexity  of  the  model,  the  concentra¬ 
tions  of  reactant  gases  in  a  fuel  cell  can  initially  be  expressed  as 
the  inlet  partial  pressure,  /?n2,i  and  Po2,\  (i-e.  as  the  feed  enters 
the  particular  gas  channel)  or  the  average  partial  pressure  in 
the  channel.  The  average  channel  parameters  would,  therefore, 
be  the  hydrogen  partial  pressure  in  the  anode  gas  channel  and 
the  average  oxygen  partial  pressure  in  the  cathode  gas  channel, 
/?n2,ch  and  /?o2,ch>  respectively.  These  ‘channel’  partial  pres¬ 
sures  would,  in  general,  require  correction,  by  appropriate  mass 
transfer  calculations  in  the  MEA,  to  the  lower  values,  and 
Pq2  ,  at  the  ‘reaction  interfaces’  where  all  the  participants  in  the 
reactions  come  together.  Typically,  the  latter  should  actually  be 
expressed  as  liquid-phase  parameters,  the  dissolved  gas  concen¬ 
trations  c^?  and  Cqo  . 

Some  modelling  approaches  include  the  mass  transfer  losses 
in  a  more  visible  way,  including  concentration  polarization  terms 
for  each  electrode  and  starting  from 

V  =  E  +  77 act, a  T  0 act, c  T  7? ohmic  "f"  he one, a  T  0conc,c  (5) 

In  general,  modelling  of  77act,a  and  7/Conc,a  requires  the  prediction 
of  c from  j7H2,ch  and  the  modelling  of  r]act,c  and  r]conc,c  requires 
the  prediction  of  Cq2  from  j7o2,ch- 

4.  Modelling  reactant  concentrations  in  the  electrolyte 


4.1.  Introduction — the  exchange  current  density 


A  general  correlation  for  exchange  current  density,  essen¬ 
tially  the  electrochemical  reaction  rate,  could  take  the  form  of  a 
conventional  heterogeneous  rate  equation 


iQ  =  A  exp 


B 

T 


where  ‘A’  will  typically  contain  one  or  more  rate  constants  as 
well  as  concentration  terms  to  represent  all  the  reactants.  This 
could  take  the  general  form 


A  =  A°C^clMo]  (7) 

The  parameter  A°  will  contain  the  kinetic  constants  and  will 
normally  be  only  a  function  of  temperature.  Cm  would  be 
the  catalyst  parameter,  representing  the  total  concentration  of 
active  sites,  and  could  appear  to  a  reaction  order,  m,  other  than 
unity.  Often,  in  heterogeneous  catalysis,  A°  and  are  lumped 
together  and  are  not  quantified  individually. 

The  concentration  of  dissolved  reactant  gas,  cgas,  at  the 
electrolyte-catalyst  interface  could  also  be  raised  to  some  power 
if  the  reaction  order  with  respect  to  this  concentration  is  other 
than  unity.  In  general,  cgas  would  be  a  function  of  the  solubility 
of  the  reactant  gas  in  the  electrolyte  and  could  depend  on  the 
solute  (e.g.  H2  or  O2),  the  partial  pressure  of  the  solute,  the  sol¬ 
vent  (i.e.  the  particular  electrolyte  and  it’s  concentration),  and 
the  temperature. 

Finally,  fe  would  be  a  general  Langmuir-Hinshelwood  term, 
common  in  heterogeneous  rate  equations,  to  quantify  the  fraction 


of  the  active  catalyst  sites  that  are  ‘vacant’  and  able  to  participate 
in  the  gas-catalyst  interaction.  The/^  term  would  normally  be 
some  function  of  the  equilibrium  (zero  net  current)  fractional 
surface  coverage  of  reaction  intermediate  and  would  normally 
be  a  function  of  cgas  and  T. 

The  prediction  and  modelling  of  cgas  is  the  focus  of  this  paper. 

4.2.  Solubilities  of  reactant  gases  in  the  electrolyte 
4.2.1.  Introduction 

Published  data  on  exchange  current  densities  and  polariza¬ 
tion  come  from  studies  using  a  variety  of  electrolytes,  dilute 
(0.1  or  0.05  M,  for  example)  H2SO4  probably  being  the  most 
common.  The  major  modelling  interest  for  PEMFC,  however, 
is  the  prediction  of  electrochemical  parameters  for  all  the  way 
from  PEMFC  designs  where  the  electrolyte  phase  in  direct 
contact  with  the  Pt  catalyst  was  essentially  ‘reactant  gas  dis¬ 
solved  in  water’  to  designs  where  the  electrolyte  phase  was 
‘reactant  gas  dissolved  in  humidified  Nation™’.  In  order  to 
make  maximum  use  of  published  work  on  electrodes,  we  must 
be  able  to  evaluate  the  previously  defined  A°  constant.  This 
requires,  according  to  Eq.  (7),  the  evaluation  of  cgas,  cu2  at  the 
electrolyte/anode  interface  and  co2  at  the  electrolyte/cathode 
interface.  Similarly,  in  order  to  properly  utilize  these  A°  val¬ 
ues  for  modelling  of  PEMFC,  we  must  be  able  to  predict 
the  values  of  ch2  and  c o2  in  the  membrane  electrolyte  at  the 
gas-membrane  interface  and  at  the  operating  conditions  of  the 
fuel  cell. 

The  solubility  of  the  gaseous  reactant  in  the  electrolyte 
(i.e.  the  Henry’s  Law  linkage  of  pu2  and  ch2  or  po2  and 
co2)  will  affect  the  prediction  of  activation  polarization  and 
the  diffusion  rate  of  dissolved  reactant  gas  through  the  elec¬ 
trolyte  will  be  a  factor  in  the  prediction  of  concentration 
polarization. 

In  a  relatively  simple,  1-D,  model,  the  average  hydrogen 
and  oxygen  partial  pressures  in  the  respective  feed  channels, 
PH2,ch,av  and  j7o2,ch,av?  are  used  to  quantify  the  feed  concentra¬ 
tions.  In  a  more  complex  2-D  or  3-D  model,  the  point  values 
of  Pn2,ch  and  po2,ch  are  used  to  quantify  the  local  reactant 
partial  pressure.  Assuming  that  concentration  polarization  is 
estimated  in  a  separate  calculation,  these  /?H2,ch,av  and  /?o2,ch,av 
or  pn2,ch  and  po2,ch  become  the  effective  reactant  partial  pres¬ 
sure,  p^n  and  j7q2,  at  the  interface  of  the  Pt  catalysts  and  the 
surrounding  aqueous  medium.  At  some  point  in  the  system,  the 
‘gas-aqueous  interface’  (subscript  i),  hydrogen  and  oxygen  gas 
go  into  solution.  The  dissolved  reactant  concentrations  at  this 
interface,  cn2,i,  and  co2,i,  can  now  be  estimated  from  the  appro¬ 
priate  Henry’  Law  expression. 

The  following  form  of  Henry’s  Law  will  be  used 


with  c  and  p  being  appropriately  subscripted  for  the  particular 
interface,  typically  subscript  ‘i’  for  the  gas-‘aqueous’  interface. 
The  ‘aqueous’  phase  could  be  the  dilute  acid  electrolyte  in  a 
single-electrode  study  or  the  polymer  membrane  electrolyte  in 
a  PEMFC  modelling  application. 
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4.2.2.  Quantification  of  ‘p’ 

A  couple  of  caveats  should  be  kept  in  mind. 

Many  published  studies  report  parameter  values  (such  as  iQ, 
for  example)  “at  atmospheric  pressure”.  Ideally  this  should  mean 
‘at  pu2  =  1  atm’  or  ‘at  po2  =  1  atm’.  In  fact  it  could  simply 
mean  than  the  data  were  collected  in  an  apparatus  open  to  the 
atmosphere  so  that  ‘atmospheric  pressure’  was  the  ambient  pres¬ 
sure  at  that  time,  probably  an  absolute  pressure  between,  say,  99 
and  103  kPa,  i.e.  an  applied  pressure  of  perhaps  1  atm  db  2%. 

A  second  concern  is  the  correction  (or  lack  of  correction)  for 
the  partial  pressure  of  water  vapour.  With  ideal  behaviour  at  the 
low  pressures  typical  of  fuel  cell  operation,  we  may  assume,  for 
pure  feeds  of  reactant  gas, 

Ph2  or  P(h  =  /’total  -  Ph2o  (9) 

and 

PHiO  =  Pnlo  (!°) 

For  the  typical  temperature  range  of  interest  to  PEMFC,  say 
0-90  °C,  p^lo  will  range  from  0.006  to  0.69  atm,  values  ranging 
from  ‘insignificant’  to  potentially  ‘very  significant’.  The  relative 
humidity,  RH,  will  typically  range  from  perhaps  40  to  80%  in  a 
lab  environment,  probably  close  to  100%  in  the  gas  space  above 
an  aqueous  electrolyte  in  a  single-electrode,  perhaps  25-100% 
in  the  anode  gas  channel  of  a  PEMFC  and  close  to  100%  in  the 
cathode  gas  channel  of  a  PEMFC.  Clearly,  therefore,  there  will 
be  situations  where  a  device  run  ‘at  atmospheric  pressure’  could 
have  a  pgas  ( pu2  or  Po2 )  as  low  as  about  0.3  atm,  drastically 
different  from  ‘  1  atm’ . 

4.2.3.  Quantification  of  H 

4.2. 3.1.  Henry's  Law  constants  for  aqueous  acid  electrolytes. 
The  most  common  electrolytes  in  contact  with  PEM  electrodes 
are  water  and  water-soaked  Nation  but  published  electrode  polar¬ 
ization  data  typically  are  for  various  hydrogen  or  oxygen  elec¬ 
trodes  where  the  Pt  catalyst  is  in  contact  with  an  aqueous  acid 
electrolyte.  Solubility  data  for  hydrogen  and  oxygen  in  a  variety 
of  acid  electrolytes  are,  therefore,  necessary,  initially  to  inter¬ 
pret  and  correlate  published  values  of  exchange  current  density 
and  finally  to  apply  such  iQ  correlations  to  the  prediction  of 
activation  polarizations  in  modelled  PEMFC  at  the  operating 
temperature  and  reactant  compositions  of  the  fuel  cell.  Recom¬ 
mended  expressions  for  Henry’s  Law  ‘constants’  for  various 
systems  are  summarized  in  the  following  sections. 

4.2. 3. 1.1.  Henry’s  Law  constants  for  hydrogen  in  aqueous 
acid  electrolytes.  For  H2  in  water:  Analysis  of  data  from  var¬ 
ious  sources  [15-19]  led  to  the  following  proposed  correlation 
for  temperatures  from  0  to  100  °C  and  pu2  up  to  10  atm: 

from0to45°C  :  //h2,h2o  =  7.9  x  106exp(— 545/T) 

X  [1  +0.000071/42],  (11) 

from45to  100°C  :  //h2,h2o  =  8-34  x  105exp(170/7’) 

x  [1  +  0.000071/7^].  (12) 


For  H2  in  dilute  H2SO4  electrolyte:  Analysis  of  25  °C  data 
from  the  International  Critical  Tables  ([1],  p.  145,  Fig.  24)  led 
to  the  recommendation  that  the  above  ‘pure-water’  value  can  be 
multiplied  by  the  following  concentration  correction  factor  for 
up  to  about  4N  acid  (i.e.  [N]  =  4)  or  about  2  M  acid  (i.e.  [M]  =  2): 

7/h2,h2so4  =  T/h2,h2o(1  +  0.1[A])  (13) 

7/h2,h2so4  =  7/h2,h2o(1  +  0.2[M])  (14) 

Lacking  any  data  showing  temperature  effects,  it  has  been 
assumed  that  the  temperature  effects  for  the  ^/dilute  H2SO4 
system  are  the  same  as  those  for  the  H2/H2O  system  and  are, 
therefore,  already  contained  in  the  correlations  recommended 
above.  Note  that  the  commonly  used  0.05  and  0.1  M  H2SO4 
electrolytes  would,  therefore,  have  Hu2  values  only  about  1-2% 
above  those  for  water  so  that  the  corresponding  ch2  values  would 
only  be  1-2%  lower  in  the  dilute  acids  than  in  water. 

For  H2  in  dilute  HCl  or  dilute  HNO3  electrolyte:  The  same 
ICT  reference  as  above  [1]  contained  data  indicating  that  Eqs. 
(13)  and  (14)  for  dilute  H2SO4  should  be  sufficiently  accurate  for 
these  dilute  acids  as  well.  For  greater  precision,  //h2,hci  could 
be  ^5%  below  //h2,h2S04  and  //h2,hno3  could  be  ^8%  below 
//h2,h2so4- 

For  H2  solubility  in  concentrated  H3PO4 :  Data  from  PAFC 
studies,  although  in  a  different  acid  and  generally  at  temperatures 
well  above  PEM  operating  temperatures,  could  be  applicable  to 
PEM  modelling.  Although  data  for  H2  have  not  yet  been  located, 
solubility  data  for  O2  in  aqueous  acids  [see  below]  could  imply 
#h2,h3po4  values  ~20%  below  7/h2,h2so4- 

4.2.3. 1.2.  Henry ’s  Law  constants  for  oxygen  in  aqueous  acid 
electrolytes.  For  O2  in  water:  Correlating  data  from  various 
sources  [15-19],  the  following  are  recommended: 

from0to45°C  :  //o2,h2o  =  1.34  x  108  exp[— 1540/T]  (15) 

from 45  to  100  °C  :  7/o2,h2o  =  5.08  x  106  exp[— 500/T] 

(16) 

For  O2  in  dilute  H2SO4  electrolyte:  Correlating  solubility 
data  at  25  °C  [1,15],  the  following  concentration  correction  fac¬ 
tors  were  obtained  for  concentrations  below  2N  (i.e  [N]  =  2)  or 
1  M  (i.e.  [M]  =  1): 

7/o2,h2so4  =  7/o2,h2o(1  +  0.125[A])  (17) 

T/o2,h2so4  =  T/o2,h2o(1  +  0.25[M])  (18) 

Lacking  any  data  showing  temperature  effects,  it  has  been 
assumed  that  the  temperature  effects  for  the  02/dilute  H2SO4 
system  are  the  same  as  those  for  the  O2/H2O  system  and  con¬ 
tained  in  Eqs.  (15)  and  (16)  correlations. 

For  O2  solubility  in  H3PO4 :  Published  data  ([1],  p.  144,  Fig. 
43)  for  solubilities  up  to  about  85  wt.%  H3PO4  have  been  utilized 
as  follows: 

For  the  dilute  region:  For  wt.%  (C)  from  0  to  60, 
7/o2,h3po4,25°c  =  0-78  x  106  exp(0.0157C)  (19) 

For  ~  100  wt.%  H3PO4,  7/o2,h3po4,25°c  =  6.9  x  106 
(extrapolated  estimate). 
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For  O2  in  dilute  HCl  or  dilute  HNO3  electrolyte:  From 
ICT  ([1],  p.  145,  Fig.  44),  these  solubilities  at  25  °C  lie  in 
between  those  of  O2  in  H2SO4  and  O2  in  dilute  H3PO4,  Z7h2,hci 
being  ~ 5 %  above  Z7o2,h2so4  and  Z7o2,hno3  being  about  10% 
above  Z7o2,h2so4-  Since  the  earlier  recommended  correlations 
for  //o2,h2so4  included  the  assumption  of  the  same  temperature 
variation  as  Z7o2,h2o>  that  assumption  is  automatically  carried 
forward  to  these  acids  as  well. 

4. 2. 3. 2.  Henry’s  Law  constants  for  Nafion™  membrane  elec¬ 
trolytes.  There  are  a  number  of  ion-exchange  membranes  which 
have  been  applied  to  the  development  of  ‘proton-exchange  mem¬ 
brane’,  PEM,  fuel  cells.  Nafion™  is  the  membrane  of  major 
interest  in  the  present  work  although  reference  will  be  made  to 
other  related  materials. 

4.2. 3. 2.1.  Henry’s  Law  constants  for  hydrogen  in  Nafion. 
Yeo  and  McBreen  [20]  are  often  cited  as  an  early  source  of 
hydrogen  solubility  data  in  Nafion. 

Ogumi  et  al.  [21],  in  their  major  report  on  oxygen  solubility 
in  Nafion  120  (see  Section  4.23.2.2  for  more  details),  include 
a  25  °C  value  of  24.50  x  10~3  M  for  the  solubility  of  hydrogen 
in  Nafion  120,  apparently  from  Yeo  and  McBreen  [20].  This 
probably  is  equivalent  to  an  /7ho  ,Nafion  value  of  about  4.08  x  104. 

Sakai  et  al.  [22]  reported  solubilities  of  H2  in  dried  Nafion 
125,  acid  H+-form  and  K+-form.  In  the  0-45  °C  region,  the  two 
sets  of  data,  when  converted  to  Henry’s  constants,  had  Arrhenius 
temperature  coefficient  terms  of  exp(— 363/7)  and  exp(— 41 1/7), 
respectively.  These  compare  reasonably  with  the  corresponding 
#h2,h2o  coefficient  of  exp(— 545/7).  From  45  to  100  °C,  there 
were  only  K+-form  data  and  the  temperature  coefficient  term  is 
exp(— 1040/7)  with  the  corresponding  term  for  Z7h2,h2o  being 
exp(+ 170/7).  The  Nafion  appears  to  behave  water-like  below 
45  °C  but  not  above  that  temperature.  At  the  typical  temperatures 
for  PEMFC  operation  (the  45-100  °C  range),  Z7h2,h2o  values 
appear  to  be  from  3^  to  6^  times  Z7Ho,  Nafion  1 25, k+ -dried  values, 
the  latter  being  given  by: 

#H2,  Nafion  125, K+ -dried  =  5.6  X  106  exp(- 1040/  T)  (20) 

The  Henry’s  Law  recommendations  by  Bernardi  and  Ver- 
brugge  [23]  were:  “Yeo  and  McBreen  [20]  report  the  solubility 
of  hydrogen  at  25  °C  in  Nafion  120.  Data  at  higher  tempera¬ 
tures  were  not  reported.  Since  the  solubility  of  hydrogen  in 
water  is  only  weakly  dependent  on  temperature  in  the  25-100  °C 
temperature  range,  the  Henry’s  constant  is  assumed  to  be  inde¬ 
pendent  of  temperature.”  (Support  for  this  last  statement  is  a 
reference  to  ‘Berger,  p.  143’  [1,2].  See  above  for  further  detail 
re  this  reference.)  Using  this  argument,  their  ZZh2, Nafion, 25  °c  of 
4.5  x  104  atm  cm3  mol-1  should  be  combined  with  the  temper¬ 
ature  dependence  for  Z7h2,h2o  given  by  Eqs.  (11)  and  (12)  to 
give,  for  the  45-100  °C  temperature  range: 

Hu2, Nafion  =  0.255  x  105  exp(170/  T)  (21) 

i.e.  only  ~3%  of  77h2,h2o 

4. 2. 3. 2. 2.  Henry’s  Law  constants  for  oxygen  in  Nafion. 
Ogumi  et  al.  [21]  appear  to  be  widely  cited  as  the  first  to  publish 
data  for  oxygen  solubility  in  Nafion.  They  used  Nafion  120,  with 


0.83  meq  g-1  of  exchange  capacity,  in  contact  with  “prehumid¬ 
ified  oxygen  ...  at  atmospheric  pressure”.  It  is  not  completely 
clear  what  the  value  of  po2  is  for  their  reported  solubility  results 
since,  for  ‘prehumified  gas’, 

Po2  =  ^totai  -  Pu2o  (22) 

and  pu2o  could  have  any  value  from  0  up  to  /?H2o,sat  f°r  the 
particular  test  temperature.  This  does  give  a  range  for  pu2o  ancf 
therefore,  a  range  for  the  calculated  values  of  77o2, Nafion-  The 
lower  values  of  the  77- values  below  are  probably  the  more  likely. 


T(°C) 

c„  x  103  (M) 

^02, Nafion  X 

10-5  (atm  cm3  mol-1) 

20 

7.20 

1.36-1.39 

30 

6.50 

1.47-1.54 

40 

5.30 

1.75-1.89 

50 

5.90 

1.49-1.69 

Their  membranes  were  stated  to  be  “pretreated  by  soaking  in 
boiling  water  in  the  Na+  form  for  30  min”.  There  was  also  the 
observation  reported  that  “oxygen  solubility  was  almost  inde¬ 
pendent  of  the  sample  and  pretreatment”. 

Ogumi  et  al.  [24],  for  both  oxygen  and  hydrogen  in  Nafion 
125,  stated  that  “solubility,  c°,  was  very  high  compared  with  the 
value  in  aqueous  solutions,  was  of  the  order  10-2  mol  dm-3,  and 
was  almost  independent  of  the  water  content  of  the  membrane.” 
Assuming  that  the  same  experimental  apparatus  and  procedure 
were  used  as  in  the  1984  work  [21],  there  is  a  small  uncertainty 
in  the  actual  value  of  po2  when  converting  the  solubility  to  a 
value  of  77o2, Nafion-  At  the  25  °C  used  for  the  following  results, 
the  actual  77  could  be  1  above  or  below  the  value  shown: 


Nafion  form 

•^02,  Nafion 

(atm  cm3  mol-1) 

Nafion  125-Na 

7.2  x  104 

Nafion  117-Na 

1.14  x  105 

Nafion  125-K 

9.03  x  104 

Nafion  117-K 

9.21  x  104 

Several  relevant  ratios,  for  25  °C,  can  also  be  calculated  from 
their  results: 


#H2,H2o/#H2,  Nafion  120  =  0.032, 

#h2,h2so4/#h2,h2o  =  1-08, 

#H2, Nafion  1 20 /#02, Nafion  120  ~  1,  and 

#o2,h2o/#h2,h2o  =  0.6. 

Sakai  et  al.  [25]  report  extensive  data  on  O2  and  H2  perme¬ 
abilities  through  Nafion  117  and  125  but  make  no  direct  mention 
of  solubility  data. 

Sakai  et  al.  [22]  reported  solubilities  of  O2  in  dried  Nafion 
125,  acid  H+-form  and  K+-form,  Nafion  117-K+(dry)  and 
Nafion  117-H+(wet).  The  four  sets  of  data,  when  converted  to 
ZZo2,Nafi0n-xx  and  put  on  an  Arrhenius  plot,  lay  in  four  essen¬ 
tially  parallel  lines,  temperature  coefficient  exp(— 970/7).  All 
the  ‘dried  membrane  data,  all  for  25-100  °C,  were  fitted  by 

tfo2,Nafion-xx  »  3.73(±25%)  x  105  exp(-970/  T)  (23) 
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while  the  few  ‘wet  membrane’  //-values,  all  taken  from  data  at 
or  below  25  °C,  were  50-80%  higher  than  values  given  by  the 
extrapolated  ‘dried  membrane’  correlation. 

Bernardi  and  Verbrugge  [26]  cite  the  solubility  data  published 
by  Ogumi  et  al.  [21]  (and  tabulated  above)  as  the  basis  of  their 
derived  Henry’s  constant  correlation: 

tfo2,Nafion  =  1-33  x  106  exp(-666/ T)  (24) 

Strictly  speaking,  this  correlation  was  based  on  data  from 
20  to  50  °C  only  and,  as  is  clear  in  the  above  tabulation,  the 
Arrhenius  behaviour  of  the  40  and  50  °C  data  is  a  bit  inconsistent 
with  the  data  below  40  °C.  Combining  this  with  observations  in 
other  publications  that  there  is  commonly  an  inflection  point  in 
the  40-50  °C  range  of  an  Arrhenius  plot  of  Nation  membrane 
characteristics  (water  content,  dissolved  gas  content,  etc.),  the 
above  Bernardi  and  Verbrugge  correlation  probably  should  not 
be  used  above  40  °C.  Nevertheless,  their  Table  3  includes  an 
80  °C  value  consistent  with  Eq.  (24)  correlation.  In  a  subsequent 
paper  [23],  Bernardi  and  Verbrugge  quote  their  1991  paper  but 
suggest  the  correlation 

In  Ko2  =  -(666/  r)  +  14.1  (25) 

which,  since  H  and  Ko2  have  identical  definitions,  is  consis¬ 
tent  with  their  1991  correlation.  However,  their  1992  Table  II 
[23]  includes  two  values  of  Henry’s  constant:  an  80  °C  value 
of  2.0  x  105  (which  is  consistent  with  the  above  Eqs.  (24)  and 
(25))  and  a  95  °C  value  of  8.1  x  106  (which  is  far  from  being 
consistent  with  Eqs.  (24)  and  (25)). 

Parthasarathy  et  al.  [27,28]  reported  extensive  data  on  water 
content  and  oxygen  solubility  in  Nation  from  25  to  80  °C  and  at 
an  oxygen  pressure  of  5  atm.  The  //o2,Nafion  values  range  from 
3.12  x  105  at  25  °Cto  1.13  x  106  at  80° C  with  a  major  change  in 
slope  on  an  Arrhenius  plot  at  8.7  x  105  and  40  °C.  They  do  com¬ 
ment,  however,  that  “the  high  solubilities  (i.e.  low  //o2,Nafion) 
at  low  temperature  (<50  °C)  do  not  persist  after  the  first  heat¬ 
ing  cycle.  After  the  polymer  has  imbibed  sufficient  water  at  the 
higher  temperature,  lowering  the  temperature  does  not  enhance 
the  oxygen  solubility.  This  indicates  that  water  imbibed  by 
the  Nation  is  permanently  associated  with  the  structure  of  the 
polymer.”  This  conclusion  suggests  the  possibility  that  many 
published  solubility  data  are  at  Nation  conditions  not  represen¬ 
tative  of  the  operating  conditions  of  a  PEM  fuel  cell  and  raises  a 
warning  flag  about  the  interpretation  and  use  of  published  solu¬ 
bility  data.  Subsequently,  Gottesfeld  and  Zawodzinski  [29]  state 
that  the  above  solubility  results,  utilizing  Pt/bulk  ionomer  rather 
than  Pt/recast  ionomer,  ‘do  not  mimic  precisely  the  interfacial 
compositions  in  a  Pt/C/recast  ionomer  composite’. 

Broka  and  Ekdunge  [30]  reported  results  from  extensive 
oxygen  and  hydrogen  permeation  and  water  uptake  studies  on 
Nation  membrane  but  did  not  report  any  gas  solubility  results. 

Gottesfeld  and  Zawodzinski  [29]  include  (in  their  Table  2, 
p.  205)  some  solubilities  which  permit  the  following  semi- 
quantitative  comparisons: 

^02 ,  recast  Nafion  ionomer  film/#O2,0.5MH2SO4  ~  0-36 
^02,  bulk  Nafion  ionomer /^02, 0.5  MH2S04  ^  0.07—0.13 


This  clearly  shows  the  large  difference  between  the  bulk  and 
recast  Nafion  ionomer  but  also  can  be  used,  in  conjunction  with 
our  earlier  correlations  for  //o2,h2so4  and  //o2,h7o>  to  suggest 
//o2, Nafion  correlations  that  are  perhaps  more  applicable  to  mod¬ 
ern  PEMFC: 

from  0  to  45  °C  :  H0l,  recast  Nafion  ~  5.43  X  107  exp[— 1540/  T] 

(26) 

from  45  to  100  °C  :  //o2, recast  Nafion  ^2.06  x  106  exp[  — 500/  T] 

(27) 

It  is  worth  noting  that  these  ‘synthesized’  //o2, Nafion  correla¬ 
tions  lie  in  the  middle  of  the  Arrhenius  plot  of  all  the  literature 
data  summarized  in  the  present  paper  and  about  a  factor  of  2^ 
above  the  often-used  1991  and  1992  correlations  of  Bernardi 
and  Verbrugge  [26,23];  that  is,  predicted  oxygen  solubilities  in 
Nafion  would  be  about  40%  of  the  Bernardi  and  Verbrugge  pre¬ 
dictions.  Such  a  40%  difference  in  co2, Nafion  would  effectively 
mean  a  40%  difference  in  the  predicted  exchange  current  densi¬ 
ties  at  the  cathode. 

Um  et  al.  [31]  quantified  the  oxygen  solubility  in  Nafion 
membrane  at  80  °C  by  an  //o2, Nafion  value  of  ~2  x  105  atm  cm3 
mol-1.  This  appears  to  be  simply  an  acceptance  of  the  above- 
mentioned  recommendations  of  Bernardi  and  Verbrugge 
[26,23]. 

Mitsushima  et  al.  [32]  reported  considerable  solubility  data 
for  four  membranes  of  different  exchange  capacities,  one  of 
them  being  Nafion  1 17  at  ^0.91  meqg-1,  at  a  po2  of  5  atm.  As 
noted  elsewhere,  there  is  the  uncertainty  as  to  whether  or  not  the 
partial  pressure  of  water  vapour  has  been  allowed  for  in  stating 
the  value  for  po2  but  this  uncertainty  is  minor  at  5  atm.  The 
authors  do  state,  however,  that  “measurements  carried  out  with 
the  feed  gas  at  saturated  relative  humidity”.  The  solubilities,  on 
conversion  to  Henry’s  constants,  are: 


r(°c) 

^02,  Nafion 

(atm  cm3  mol-1) 

20 

7.24  x  105 

30 

7.8  x  105 

50 

1.54  x  106 

70 

5.55  x  106 

A  number  of  PEM  modelling  papers  by  Broka  (Dannenberg) 
et  al.  [30,33,34]  include  considerable  data  for  O2  permeability 
in  Nafion,  “based  on  solubility  measurements”  but  only  say  “the 
solubility  of  the  oxygen  is  higher  in  PTFE  than  in  water”.  They 
refer  to  several  earlier  publications  [22,24,25],  all  discussed 
above. 

Comments  and  conclusions  re  //o2, Nafion-  If  all  the  above- 
cited  results  are  put  on  an  Arrhenius  plot,  a  huge  spread  in 
the  published  numerical  values  of  //o2,  Nafion  is  apparent,  rang¬ 
ing  from  about  an  order  of  magnitude  at  low  temperatures  to 
approaching  two  orders  of  magnitude  at  100  °C.  There  are  prob¬ 
ably  several  contributing  factors  to  this  apparent  major  obstacle 
to  the  accurate  prediction  of  co2,  Nafion  in  modelling  of  PEMFC. 

Gottesfeld  and  Zawodzinski,  in  their  major  publication  on  all 
aspects  of  PEMFC  [29],  dwelt  at  length  on  the  importance  of 
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using  data  for  the  appropriate  form  of  Nation,  some  properties 
of  bulk  ionomer  being  significantly  different  than  those  for  the 
recast  ionomer  commonly  used  in  PEM  ME  As. 

The  appropriate  water  content  in  the  membrane  is  also  very 
important  as  gas  solubilities  are  certainly  dependent  on  the  actual 
‘imbibed  water/Nafion  matrix  system’  used  to  measure  solubil¬ 
ities  on  the  one  hand  or  to  model  particular  PEMFC  on  the  other 
hand.  The  recent  work  by  Fimrite  et  al.  [35],  although  containing 
no  solubility  data,  gives  a  good  overview  of  a  modelling  frame¬ 
work  for  the  evaluation  of  physical  properties  in  water/Nafion 
systems. 

The  solubility  data  in  the  literature  could  also  be  subject  to 
a  variation  based  on  the  experimental  details,  often  not  consid¬ 
ered  important  enough  to  report  completely.  The  preparation  of 
the  membrane  prior  to  each  run,  the  humidification  procedure 
followed,  the  sequence  of  the  experimental  runs  (for  example, 
was  a  single  membrane  sample  tested  at  progressively  increas¬ 
ing  temperatures  or  progressively  decreasing  temperatures  so 
that  the  water  contents  suffered  from  a  ‘memory’  of  the  previ¬ 
ous  water  content)  and  the  duration  of  each  run  are  all  possibly 
significant.  We  have,  for  example,  observed  PEM  cells  taking 
days  rather  than  hours  to  reach  equilibrium  (as  evidenced  by  a 
completely  constant  voltage),  a  feature  that  we  attributed  to  slow 
equilibration  of  the  water  content  of  the  membrane. 

5.  Summary  and  conclusions 

In  general,  prediction  of  activation  and  concentration  polar¬ 
ization  in  a  PEMFC  model,  should  pay  more  attention  to  the 
proper  evaluation  of  pgas  and  cgas. 

Expressions  such  as  Eqs.  (6)  and  (7)  must  be  derived  to 
provide  iQ,  essentially  the  kinetic  rate  equation  for  the  rate¬ 
determining  step(s).  This  expression  will  require  surface  con¬ 
centrations  that  will  be  modelled  using  expressions  such  as  Eqs. 
(8)-(10)  and  (1 1)— (27). 

A  major  modelling  challenge  is  the  selection  of  the  appropri¬ 
ate  Henry’s  Law  correlation  for  the  ME  A  being  modelled,  the 
form  and  the  water  content  of  the  Nafion  at  the  gas/Pt/Nafion 
interface  having  a  large  influence  on  the  resultant  concentration 
of  the  reactant  gas  at  the  catalyst  surface. 
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